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ABSTRACT 

We study the arc statistics of gravitational lensing generated by dark matter 
halos in order to probe their density profile. We characterize the halo profile by 
two parameters, the inner slope of the central cusp a, and the median amplitude 
of the concentration parameter, c norm , for a halo of mass 10 1A h~ 1 M Q at z — 0, 
and compute the numbers of tangential and radial arcs produced by gravitational 
lensing of galaxy clusters. We find that the number of arcs divided by the number 
of halos is a good statistic which is sensitive to both c norm and a with very weak 
dependence on the cosmological parameters. If the arc samples with well-defined 
selection criteria for the clusters become available, one can strongly constrain 
both c n orm and a. While our tentative comparison with the existing observational 
data indicates that the inner density profile of dark halos is indeed as steep as 
predicted by recent simulations (a ~ 1.5), the homogeneous samples of tangential 
and radial arcs are required for more quantitative discussions. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general 
— gravitational lensing 

1. Introduction 

Dark matter halos play a central role in the standard picture of the cosmological struc- 
ture formation as plausible sites hosting a variety of astronomical objects such as galaxies 
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and clusters of galaxies. On the basis of a series of systematic cosmological N-body simula- 
tions, Navarro, Frenk, & White (1996, 1997, hereafter NFW) found that the density profile 
obeys the "universal" form p{r) oc r _1 (r + r s )~ 2 irrespective of the underlying cosmological 
parameters, the shape of the primordial fluctuation spectrum and the formation histories. 
More recent high-resolution simulations indeed confirmed the existence of the central cusp 
but suggested the even steeper inner slope; p{r) oc r _L5 rather than oc r _1 (Moore et al. 
1999; Fukushige & Makino 2001). Nevertheless the universality of the profiles in numerical 
simulations is fairly established except for the possible weak dependence on the halo mass 
(Jing & Suto 2000) and also for some scatter around the mean (Jing 2000). 

The above indications from simulations, however, do not seem to be supported by ei- 
ther simple theoretical considerations or available observations. Plausible theoretical models 
rather predict that the inner slope of the halo profile should depend on the primordial fluc- 
tuation spectrum (Hoffman & Shaham 1985; Syer & White 1998) and also on the merging 
history (Nusser & Sheth 1999; see also Lokas & Hoffman 2000). Detailed analyses of the 
X-ray surface brightness of clusters of galaxies (Wu & Xue 2000; Wu & Chiueh 2001) are 
inconsistent with the scaling of the halo concentration against the halo mass predicted by 
simulations. Furthermore both rotation curves of the low surface brightness galaxies (Salucci 
& Burkert 2000; de Blok et al. 2001) and the inner region of the cluster CL0024-1654 recon- 
structed from gravitational lensing images (Tyson, Kochanski, & Dell'Antonio 1998) indicate 
the flat core instead of the central cusp, although some controversy about each claim still 
remains (van den Bosch et al. 2000; Broadhurst et al. 2000; Shapiro & Iliev 2000). This 
conflict has motivated wild proposals, including an idea that dark matter is self-interacting 
(Spergel & Steinhardt 2000). 

Since the current situation concerning the numerical, theoretical and observational in- 
dications for dark matter halo profiles is somewhat puzzling, it is important to develop 
another independent methodology to probe the profiles. For this purpose, we focus on the 
arc statistics of the gravitational lensing in the present paper. The major advantages of 
this methodology include, (i) the gravitational lensing offers us the direct route to the mass 
distribution of the dark halo without additional assumptions, for instance, on the physical 
conditions of gas and stars, (ii) the gravitational arcs are produced mainly due to galaxy clus- 
ters which have an empirically good one-to-one correspondence with dark halos, in marked 
contrast with the case of the multiple QSO images due to galaxies (Nakamura & Suto 1997; 
Li & Ostriker 2001; Wyithe, Turner, & Spergel 2001; Keeton & Madau 2001; Takahashi & 
Chiba 2001), and (iii) observational confrontation for an individual object may suffer from 
the specific selection function and the scatter from the mean profile, and thus the statistical 
average over the cosmological volume is important. 
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Several authors have already examined the effect of the inner profile of dark halo on giant 
luminous arcs (Wu & Hammer 1993; Miralda-Escude 1993a,b; Hamana & Futamase 1997; 
Hattori, Watanabe, & Yamashita 1997; Williams, Navarro, & Bartelmann 1999; Molikawa 
et al. 1999; Meneghetti et al. 2001). In particular Bartelmann et al. (1998) suggested a 
strong dependence of arc statistics on cosmological parameters because of the different core 
structure of dark halos for different cosmological model. Keeton & Madau (2001) found 
that the number of predicted lenses is strongly correlated with core mass fraction, which 
results in a strong degeneracy between the inner slope of the central cusp and the dark 
matter concentration. Molikawa & Hattori (2001) pointed out that the number ratio of 
tangential and radial arcs produced by a given cluster is tightly correlated with the inner 
slope; on the basis of eleven tangential and three radial arcs for six clusters, they conclude 
that the central cusp oc r~ L5 is indeed favored. Although the existing samples of clusters are 
somewhat heterogeneous and do not satisfy well-defined selection criteria, this indicates that 
the arc statistics become useful probes of the core structure of the dark halos. Therefore we 
present a first systematic study of the effects of the dark halo profiles on the gravitational 
tangential and radial arc statistics. In particular we take proper account of the magnification 
bias, the finite size of the source galaxies, and the luminosity distribution and evolution of 
source galaxies. 

The rest of the paper is organized as follows. Section 2 briefly summarizes the main 
properties of the generalized NFW profile, and §3 presents an analytic formalism of the 
number count of arcs. Our predictions for the arc statistics are shown in §4. Finally we 
summarize the conclusions and discuss further implications in §5. 



2. Description of the Density Profiles of Dark Matter Halos 

2.1. Generalized NFW Profile 

Throughout the paper, we adopt a generalized NFW profile for dark matter halos of a 
form (Jing & Suto 2000): 

Pcrit^c 

(r/r s ) a (1 +r/r s y 

where r s is a scale radius and 5 C is a characteristic density. The profile with a = 1 corresponds 
to the one NFW proposed, and that with a — 1.5 agree with the inner profile claimed by 
Moore et al. (1999) and Fukushige & Makino (2001). The shape of halos is also characterized 
by the concentration parameter, which is defined as the ratio of the size of the halo to the 
scaled radius r s . Originally NFW used r 2 oo, the radius of halo where the mean inner density 
reaches 200 times the critical density of the universe. Rather we follow Bullock et al. (2001) 
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and adopt the definition: 

The virial radius, r vir , in the above expression is defined through the overdensity A vir at the 
virialization epoch 

M vir = —A vir p(z vir )r 3 viT , (3) 

where p(z vir ) denotes the mean density of the universe at virialization. Once the density pa- 
rameter fl and the cosmological constant A are specified, the value of A vir can be computed 
using the nonlinear spherical collapse model. We use the following formulae: 

{^ 2 (coshr] vir - l) 3 ^ 1 Ao = ^ 

(sinh?7 v i r - 7] vir ) 2 " (4) 

18tt 2 (1 + 0.4093< r 9052 ) (fi < 1, A = 1 - fi ), 

where ?7 v i r = cosh _1 (2/r2 vir — 1), u vir = 1/Q v i r — 1, and the density parameter at virialization 
is 

O = ^0(l+^vir) 3 

V " Q (l + ^vir) 3 + (1 " fi " A )(l + ^vir) 2 + A ' 1 ^ 

The approximation to A vir in equation (4) for Ao = 1 — is obtained by Kitayama & Suto 
(1996). 

Equations (1) and (3) imply that the characteristic density 5 C is related to the concen- 
tration parameter c vir as 

A • Q ■ r 3 

3 A(c v i r ) 



where A(c v i r ) is 

with 2 F\ (a, 6; c; x) being the hypergeometric function (e.g., Keeton & Madau 2001) 



c 3-a 

-4(c v ir) = „ vir 2-F1 (3 - a, 3 - a; 4 - a; -c vir ) , (7) 
3 — a: 



2.2. Concentration Parameter 

Navarro, Frenk, & White (1997) and Bullock et al. (2001) have extensively examined 
the cosmological model dependence and redshift evolution of the concentration parameter 
from N-body simulations, adopting the profile (1) with a = 1. Since we consider models 
with a 7^ 1 as well, we have to generalize their results. For this purpose, we follow Keeton & 
Madau (2001). They first define the radius r_ 2 at which the logarithmic slope of the density 
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profile is —2, i.e., d\np/dlnr = —2. For the profile of equation (1), r_ 2 = (2 — a)r s , and 
thus the corresponding concentration parameter reduces to 

r ■ 1 

' vir / \ 

C_2 = = 7: C vir- (8) 

r_ 2 2 — a 

Then they point out the importance of the scatter of the concentration parameter (Jing 
2000; Bullock et al. 2001) on the lensing statistics, and model the probability distribution 
function of c_ 2 as a log-normal function: 

(lnC_2 - lnC_ 2 ,mcdian) 5 



p(c_ 2 )dc_ 2 = —j= — exp 

V2"7RT C 



2al 



<ilnc_2, (9) 



with a c = 0.18. Jing (2000) reported that this dispersion is fairly insensitive to the cos- 
mological model parameters. Finally we introduce the scaling of c_2, m cdian according to the 
simulations by Bullock et al. (2001): 

1 / M vir \-°- 13 

C-2,median = C norm — { 1Q U h ™ Mq ) > ( 10 ) 

where h denotes the Hubble constant in units of 100 km/s/Mpc. Bullock et al. (2001) 
estimate c norm ~ 8 from their simulations with (f2o,Ao) = (0.3,0.7). In the statistical 
analyses presented below, we parameterize the halo profiles by the amplitude c norm and 
the inner slope a, assuming that the above description is applicable equally well to the three 
cosmological models that we consider. 



3. Tangential and Radial Arc Statistics 

3.1. Lens Equations 

We denote the image position in the lens plane by £ and the source position in the 
source plane by ff. For the spherical symmetric profile (1), the lens equation (e.g., Schneider, 
Ehlers, & Falco 1992) reduces to 

y = x-bf(x), (11) 

— * 

where x = |£|/r B , y = \rf\DoL/r s Dos, and Dol and Dos denote the angular diameter dis- 
tances from the observer to the lens and the source planes, respectively. The factors b and 
f(x) are related to the dark halo profile as follows: 

b = 4Pcrife , (12) 

^crit 

f(x) = - I dz I dx' - — (13) 



X 



Jo 
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where S crit is the critical surface mass density: 



_ c 2 Dps (u) 



3.2. Number of Arcs per Halo 

We consider the distortion of images of source galaxies due to the spherical halo lensing, 
neglecting the intrinsic ellipticity of those galaxies. In this case the tangential and radial 
stretching factors of a source image at the source position x with respect to the center of 
the lens halo are simply given by fJ>t(x) = (y/x)^ 1 , and fr(x) = (dy/dx)^ 1 because of the 
spherical symmetry of halos. In terms of these, we define the tangential and radial arcs as 
those satisfying 

> e t h (tangential arc), (15) 

> e t h (radial arc). (16) 

We adopt e t h = 4 for the threshold of the length-to-width ratio. This value is different from 
the canonical threshold for giant luminous arcs, e th = 10 adopted by Wu & Hammer (1993), 
because most radial arcs presently observed have a rather small length-to-width ratio. 

The projected areas around a given halo of mass M at z L satisfying equations (15) and 
(16) yield the cross sections cr(M, z L , z s ) for tangential and radial arcs of circular galaxies 
located at zs, respectively. Actually the above definition ignores the finite size of the source 
galaxies. If the size of the area on the source plane, Ai], satisfying the conditions (15) or (16) 
is smaller than that of source, however, observable arcs are not produced (e.g., Schneider et 
al. 1992). Since the smallest galaxy size that is observed as an arc (Hattori et al. 1997) 
roughly corresponds to i] crit = lh~ l kpc, we set a(M, z L , z s ) = for Ai] < i] crit . 

Figure 1 shows the lensing cross sections for tangential and radial arcs as a function of 
source redshift zs- We present the cases that the mass of the lens halo is M = 10 15 /i _1 M Q 
and M = 5 x 10 15 h~ 1 M & , because radial arcs are not formed for M ~ 1O 14 /i _1 M unless 
a or c norm is unrealistically large. The redshift of the lens halo is fixed to z^ = 0.2. These 
plots show that the lensing cross section significantly increases as a and/or c nor m become 
larger. Furthermore the cross sections for tangential and radial arcs depend on these two 
parameters rather differently. In turn, separate consideration of tangential and radial arcs 
yields useful constraints on the core structure of dark halos. 



T(x) 



H r {x) 



R(x) 



H r {x) 
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Once the relevant cross section is computed, one can calculate the number of arcs 
produced by a halo of mass M at redshift z L : 

/•zs.max C( U roo 

N(M,z L ) = dz s a(M,z L ,z s ) — (l + z s ) 3 dLn g (L,z s ), (17) 

where a(M, z^, zs) is the tangential or radial cross section in the source plane, cdt/dzs 
denotes the proper differential distance at z$: 

cdt c 1 

dz~s ~ H~o (1 + z s ) y/n (l + z s ) 3 + (l-n - A )(l + ^s) 2 + A ' 

and n g (L,zs) denotes the luminosity function of source galaxies. 

We incorporate the redshift evolution of the luminosity function of galaxies for z < 1 
adopting the empirical fit by Broadhurst, Ellis, & Shanks (1988): 



log <f)(L, z) = log 0(L, 0) + (O.lz + 0.2z 2 ) log 



0(-^max,O) 



(19) 



where L max is the bright-end luminosity corresponding to M max = — 22. + 5 log h. Therefore 
we have to restrict our consideration for source galaxies up to £s,max = 1. 

For the local luminosity function 0(L, 0), we use the Schechter function normalized to 
the two degree field (2dF) galaxy redshift survey: 

0(L,O)rfL = 0*(Ay exp ^_^^ (20) 

where 0* = 0.0169/i" 3 Mpc 3 , a = -1.28, and M* = -19.73+5 log h (Folkes et al 1999). While 
those values are derived assuming Q = 1, we compute the values in different cosmological 
models by applying an appropriate scaling so that the observed galaxy number counts versus 
their flux is unchanged. 

To evaluate equation (17), we also need the lower bound of the luminosity, L m[n , which 
depends on the magnification of arcs and the limiting magnitude of the sample. First, the 
magnification of arcs becomes 

H(x) = \fi t (x)nr(x)\ = T(x) {/i r (x)} 2 = R(x) {/i t (+)} 2 , (21) 

where x denotes the position of arcs in the lens plane. The quantity /x formally diverges in 
the case of the point source. In practice, however, /i(x) saturates at about the value which 
corresponds to the value of the position deviating from the critical curve by the source size 
(e.g., Schneider et al. 1992). Therefore we assume that all arcs are magnified by the factor 
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of e t h {f^r{xt)} 2 for tangential arc and e t h {/-^(^r)} 2 for radial arcs, where x t and x v is the 
positions of tangential and radial critical curves. Second, we use m-Q < 26.5 as the B-band 
limiting magnitude of the arc observation. The apparent magnitude can be translated to the 
luminosity if we employ the K-correction in B-band: 

K(z) = -0.05 + 2.352 + 2.55^ 2 - 4.89z 3 + 1.85z 4 (22) 

for spiral galaxies (King & Ellis 1985). Taking both effects into account, L min becomes 

T . i n-0.4(m B ,max-m*) 



e t h{/ir(^t)} 



2 



(23) 



m* = M* + 5 log ( ) + K(z s ), (24) 



in the case of tangential arcs, where D\ um is the luminosity distance. 



3.3. Total Number of Arcs 



Finally we calculate the total number of arcs by integrating equation (17) as: 



tot 



/" z L,max r 

/ dz L / 

"1,111111 "A 



M min (z L ) 



cdt 

dM N(M, z L ) n PS (M, z L )(l + z L )HnD 2 OL —, (25) 



where nps is the comoving number density of halos. We use Press-Schechter mass function 
(Press & Schechter 1974): 



n PS (M,z) 



2p(z = 0)5 (z 

7T 



M 



M 



da 



M 



dM 



cxp 



2< 



(26) 



where % is the rms of linear density fluctuation on mass scale M at z = and 5 (z) is the 
critical linear density contrast given by 



3 (12tt) 2/3 
20 D(z) ' 



(27) 



with D(z) being the linear growth rate normalized to unity at z — 0. We consider two 
selection functions, i.e., the minimum mass of integration M^^zi) for definiteness; the 
first adopts the constant minimum mass independent of z^, and the other corresponds to 
the X-ray survey with the surface brightness flux limit of S min . Throughout the paper, we 
use 0.5 — 2.0 keV band for the flux. Assuming the conventional one-to-one correspondence 
between dark halos as X-ray clusters, one can relate M min (zL) with S m i n and we use the 
relation shown in Suto et al. (2000). 
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4. Results 

4.1. Number Counts of Arcs 

In the specific examples presented below, we consider three representative cosmolog- 
ical models dominated by cold dark matter (CDM); Lambda CDM (LCDM), Standard 
CDM (SCDM), and Open CDM (OCDM) with (fi , A , h, a 8 ) = (0.3, 0.7, 0.7, 1.04), 
(1.0, 0.0, 0.5, 0.56), and (0.45, 0.0, 0.7, 0.83), respectively. The amplitude of the mass 
fluctuation, <t 8 , smoothed over the top-hat radius of 8/i -1 Mpc, is normalized so as to repro- 
duce the X-ray luminosity and temperature functions of clusters (Kitayama & Suto 1997). 

Figure 2 plots the number of halos per steradian (Top panels), tangential arcs (Middle 
panels), and radial arcs (Bottom panels) between z\,— and zl+Azl/2 with A^l = 0.05, 
for a = 1.5, c norm = 8, and z$ < 1. The triangles, open squares, and filled circles indicate 
the results for LCDM, SCDM and OCDM. The numbers of arcs in the middle and bottom 
panels are divided by the number of halos plotted in the top panels. The left and right 
panels correspond to the X-ray flux-limited (Sii m = 10~ 13 erg/s/cm 2 ) and the mass-limited 
(M min = 5 x 10 14 h~ 1 M Q ) samples. For a given mass and profile of a halo, these plots indicate 
the range of z^ which mostly contributes to the formation of arcs. Both tangential and radial 
arcs are efficiently formed around z L ~ 0.2 and z L ~ 0.1 for flux-limited and mass-limited 
samples, respectively. 

Note the strong dependence on cosmological parameters. The number of halos per 
steradian as a function of z L depends on both the volume of the universe up to z L and the 
evolution of the mass function. The former mainly explains why the halos are most abundant 
in LCDM at higher redshifts, while the latter accounts for earlier declining of halo numbers 
in SCDM (top panels in Fig. 2). Although this behavior is already well-known in the study 
of cluster abundance (Kitayama, Sasaki & Suto 1998), we emphasize that the number of arcs 
per halo also increases with the presence of the cosmological constant (middle and bottom 
panels in Fig. 2), as is pointed out by Wu & Mao (1996). This directly comes from the 
dependence of the critical surface mass density (eq. [14]). As plotted in Figure 3, the value 
of S crit is smallest in LCDM for given z$ and zl, i.e., the lensing probability is largest for a 
given halo profile. At z < 0.1, SCDM produces more abundant arcs than LCDM because 
of the higher density of halos for the same halo mass, i.e., A vir p(z v i r ) becomes the largest in 
SCDM cosmology (see eq. [4]). 

Turn next to the dependence of the halo profiles (c norm and a) on the arc statistics. 
Fi gure 4 displays the contour of arc statistics for the flux-limited sample with S\\ m — 
10~ 13 erg/s/cm 2 ; LCDM (Top panels), SCDM (Middle panels), and OCDM (Bottom pan- 
els). The left and center panels indicate the number of tangential and radial arcs per halos, 
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and the right panels plot the number ratio of radial to tangential arcs. Here we integrate 
equation (25) over the lens redshift of 0.1 < z L < 0.4 when arcs are efficiently formed as we 
discussed above. Clearly both a and c norm significantly influence the arc statistics, an order 
of magnitude more than the cosmological parameters. 

If we use the same statistics: 

W = iVtot ' rad (28) 

-^tot.tan 

proposed by Molikawa & Hattori (2001) originally for a single cluster, we confirm that the 
cosmological model dependence is extremely weak even after the statistical average over 
the redshift. In particular, the number ratio W is not so sensitive to c norm and basically a 
powerful indicator of the inner slope of the dark halo profile. The number of tangential arcs 
per halos, on the other hand, is more sensitive to c nor m (Left panels). Thus combining both 
the tangential and radial arc statistics, we can constrain both a and c norm simultaneously. 

Next consider the effects of changing the selection criterion of both halos and sources. 
Figure 5 shows how the selection criterion of lensing halos alters the prediction of arc statistics 
in the case of LCDM; from top to bottom, S\ im = 10~ 13 erg/s/cm 2 , S\ im = 10 _12 erg/s/cm 2 , 
M min = W u h~ 1 M Q , and M min = 10 15 /i _1 M o . While the different selection criterion yields a 
factor of 500 difference in the number of halos, the arc statistics per halo are fairly robust, 
and remain the powerful discriminator of the halo profile. Figure 6 displays the difference of 
the prediction by altering the threshold of the length-to-width ratio, for X-ray flux-limited 
samples, S\ im = 10~ 13 erg/s/cm 2 , in LCDM cosmology. As shown in these plots, changing the 
threshold mainly changes the number of radial arcs and consequently changes the number 
ratio W drastically. This strong threshold dependence originates from the finite source size 
effect and is not incorporated properly in Molikawa & Hattori (2001). Therefore it is clear 
that the finite size of source galaxies must be considered even in the case of the number ratio 
W. 



4.2. Uncertainties of the Predictions 

The results presented above are based on a variety of model assumptions, and we would 
like to examine the extent to which they affect the conclusions. More specifically, we focus 
on the mass function for dark halos, the size of source galaxies, and the evolution of the 
luminosity function of source galaxies. 

While the Press-Schechter mass function is widely used in various cosmological predic- 
tions, recent numerical simulations (e.g., Jenkins et al. 2001) suggest that it underpredicts 
the massive halos while overpredicts the less massive halos. Sheth & Tormen (1999) proposes 
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an empirical correction for the effect as 



n ST (M,z) = A 



2p 



2ap{z = 0) S (z) 




7T 



M 



a 



M 



da 



M 



dM 



exp 



2< 



(29) 



where a = 0.707, p = 0.3, and A = 0.322. Also we consider the cases of the twice larger 
threshold for i] crit = 2h~ x kpc and no-evolution luminosity function, i.e., (f>(L,z) = 0). 
The results separately employing the above change are plotted in Figure 7. Among them 
only the non-evolution model rather changes the total number of arcs, but this may be too 
extreme. More importantly, the ratio of the tangential and radial arcs, W, still remains un- 
changed in practice even if the luminosity evolution is neglected. The twice larger threshold 
model mainly changes the number of radial arcs, and as a result this model also changes the 
number ratio W. Thus we conclude that the arc statistics that we presented above are not 
so affected by the uncertainties of the models, and are fairly robust discriminator of the halo 
profiles if the finite size of source galaxies is correctly taken into account. 



4.3. Tentative Comparison with Observations 

While there is no homogeneous sample for the arc survey available yet that satisfies our 
selection criteria, it is tempting to make a comparison with the existing data. Luppino et al. 
(1999) present the results of imaging survey for gravitational lensing in a sample of 38 X-ray 
selected clusters of galaxies. From these clusters we choose a subsample of 13 clusters which 
satisfy the condition 0.1 < z < 0.4 and ,5(0.5 - 2.0keV) > S hm = 10~ 12 erg/s/cm 2 . Within 
these 13 clusters, 15 tangential and 2 radial arcs with e t h = 4 are reported. We attempt to 
draw cosmological implications by comparing these observational values with our theoretical 
predictions. 

Figure 8 shows our tentative comparison with observations, assuming the LCDM cos- 
mology and neglecting possible systematic errors (e.g., intrinsic ellipticies of source galaxies, 
non- sphericities in lensing halos, etc.). This result implies that dark matter halos should 
have steep inner profiles as those predicted by simulations, a ~ 1.5, but rather smaller con- 
centration, Cnorm ~ 4. While it is premature to draw definite conclusions from the present 
comparison, this analysis clearly exhibits the extent to which our current methodology puts 
useful constraints on a and c nor m separately. Note again that the number of tangential arcs 
has a strong degeneracy between a and c nor m- Thus it is important to combine the number 
of tangential arcs with the number ratio W which is mainly sensitive to a. Although the 
present example does not show the strong constraint on the concentration parameter c norm , 
a more severe constraint within ~ 10% accuracy at la level will be obtained if one uses 
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the cluster samples of iV ha i ~ 100. This would be achieved irrespective of the cosmological 
models. 



5. Discussion and Conclusions 

In this paper, we study the arc statistics of gravitational lensing produced by dark 
matter halos in order to probe their density profile. Adopting the generalized NFW profile 
(1), we describe a statistical method to predict the numbers of tangential and radial arcs 
as a function of the inner slope a and the concentration parameter c norm . We incorporate 
several realistic effects including the magnification bias, the finite size of the source galaxies, 
and the luminosity distribution and evolution of source galaxies. We find that the numbers 
of arcs sensitively depend on the values of a and c nor m- In addition, the numbers of arcs, 
if divided by the corresponding number of halos, are almost insensitive to the underlying 
cosmological parameters. Therefore they prove to be a powerful discriminator of a family of 
halo density profiles suggested by recent numerical simulations. 

Molikawa & Hattori (2001) proposed to use the number ratio of tangential and radial 
arcs W to probe the density profile. We confirm that the ratio remains a useful statistical 
measure even after taking account of the average over the cosmological mass function, the 
redshift evolution, mass-dependence and the probability distribution of the concentration 
parameter. We also show that the effects of finite source size is important even in the case 
of the number ratio W. On the other hand, W is mainly sensitive to the inner slope a, and 
we argue that the complementary information on c nor m can be obtained by combining the 
number of tangential arcs per halo. 

Our major conclusion that the cosmological model dependence of arc statistics is much 
weaker than the profile parameters seems inconsistent with the claim by Bartelmann et 
al. (1998). We note, however, that they use a different value of concentration parameter 
for different cosmological models. Thus we suspect that the claimed cosmological model 
dependence actually reflects the sensitivity to the concentration parameter that we discussed 
at length. 

The preliminary comparison with observations suggests that dark matter halos should 
have steep inner profiles. This comparison, however, is inconsistent with the mass pro- 
file of CL0024-1654 reconstructed from gravitational lensing images (Tyson, Kochanski, & 
Dell' Antonio 1998) and with the rotation curves of low-surface brightness galaxies which 
indicate a flat core (Salucci & Burkert 2000; de Blok et al. 2001). Therefore it is still pre- 
mature to draw any strong conclusions, e.g., self-interacting dark matter model (Spergel & 
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Steinhardt 2000), at this point, and it is important to put many constraints on the halo 
density profile from separate and independent analyses. 

Definitely we have to improve the present methodology by taking account of more 
realistic effects. Firstly, Bartelmann, Steinmetz, & Weiss (1995) pointed out that the number 
of arcs becomes significantly larger if the intrinsic ellipticity distribution of source galaxies 
is taken into account. This will increase both the numbers of tangential and radial arcs. 
Thus these quantitative estimate is important. Secondly, we neglect effects of the galaxies 
in clusters. In particular, the central cD galaxies move the radial arc closer to the center 
(Miralda-Escude 1995), and will affect especially the number of radial arcs. On the other 
hand, the effect of cluster galaxies seems to be small and enhances the number of arcs only 
< 15% (Flores, Mailer, & Primack 2000; Meneghetti et al. 2000), though they may affect the 
number of tangential and radial arcs differently. Finally, deviations from spherical symmetry 
of lens halos change the number of arcs (Bartelmann 1995; Molikawa et al. 1999), which will 
increase the number of tangential arcs and decrease the number of radial arcs (Molikawa & 
Hattori 2001). Therefore this effect is certainly important in discussing the number ratio W 
as well as the total number of arcs. We plan to incorporate these effects in a systematic and 
quantitative fashion which will be reported elsewhere. 

We thank Takeshi Chiba, Takashi Hamana, and Ryuichi Takahashi for useful discussions 
and Xiang-Ping Wu for his instructive comments. We also thank an anonymous referee for 
many useful comments which improved the earlier manuscript. This research was supported 
in part by the Grant-in-Aid by the Ministry of Education, Science, Sports and Culture of 
Japan (07CE2002) to RESCEU. 
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Fig. 1. — The lensing cross sections for tangential (left panels) and radial (right panels) arcs 
against the source redshift z$ for an f2 =0.3 and Ao=0.7 model. The redshift of the lensing 
halo is fixed as = 0.2. Different lines correspond to the different sets of a and c norm . Top 
and bottom panels indicate the results for the halo mass of 1O 15 /i -1 M and 5 x 10 15 /i _1 M Q , 
respectively. 
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Fig. 2. — Redshift distribution of the numbers of halos per steradian (top panels), tangential 
arcs (middle panels), and radial arcs (bottom panels). The numbers of tangential and radial 
arcs in the middle and bottom panels are divided by the number of halos plotted in the top 
panels. We adopt a = 1.5, c norm = 8, and z s < 1, and plot those numbers between z^ — Az L /2 
and zl + with Az L = 0.05. Open triangles, open squares, and filled circles indicate 

the results for LCDM, SCDM, and OCDM. The left and right panels correspond to the 
X-ray flux-limited (Sii m = 10~ 13 erg/s/cm 2 ) and the mass-limited (M min = 5 x 10 14 /i _1 M Q ) 
samples. 
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Fig. 3. — Contour of the critical surface density S crit on the zs — z\, plane for three cos- 
mological models; LCDM (top panel), SCDM (middle panel), and OCDM (bottom panel). 
The contour levels for each line are S crit = 4 (solid), 5 (long- dashed), 6 (short- dashed), 7 
(dash-dotted) and 8 (dotted), in units of 10 15 /iM Q Mpc~ 2 . 
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Fig. 4. — Predicted numbers per halo for tangential (left panels) and radial (center panels) 
arcs per halos, and their ratio W = iV to t,rad/^Vtot,tan (right panels) on a and c norm plane. We 
consider X-ray flux-limited samples with Su m = 10~ 13 erg/s/cm 2 in three representative cos- 
mological models; LCDM (top panels), SCDM (middle panels), and OCDM (bottom panels). 
The number of halo iV halo is also shown for reference. The levels of contour are 1 (thick 
solid), 10 _1 (thin solid), 10~ 2 (dash-dotted) and 10~ 3 (dotted). 
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Fig. 5. — Effect of sample selection functions on the arc statistics in the LCDM model. From 
top to bottom panels, X-ray flux-limited samples with S\ im = lCT 13 erg/s/cm 2 and Sn m = 
10~ 12 erg/s/cm 2 , and mass-limited samples with M min = iO 14 /i _1 M and M min = 1O 15 /i _1 M . 
The contour levels are the same as Figure 4. 
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Fig. 6. — Effect of the different threshold of the length-to-width ratio on the arc statistics 
for X-ray flux-limited samples with Su m = 10 _13 erg/s/cm 2 in the LCDM model; e t h = 4 (top 
panels), e t h = 7 (middle panels), and e t h = 10 (bottom panels). The contour levels are the 
same as Figure 4. 
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Fig. 7. — Sensitivity of our predictions on the adopted model parameters for X-ray flux- 
limited samples with Sn m = 10~ 13 erg/s/cm 2 in the LCDM model. First row displays our 
fiducial model (same as Figure 4). Second row adopts ust (eq. [29]) instead of n PS (eq. [26]) 
for the mass function of halos. Third row adopts 77 crit = lYr 1 kpc instead of r] crit = lh^ 1 
kpc for the cutoff size of source galaxies. Fourth row assumes no evolution in the luminosity 
function. 
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Fig. 8. — Tentative constraints on a and c norm from a sample of 13 clusters with 0.1 < z < 0.4 
and S'lim = 10~ 12 erg/s/cm~ 2 . The LCDM model is assumed. Dashed lines represent the 
relation of a and c norm which reproduce the observed number for tangential arcs and the 
number ratio of radial to tangential arcs. Dark shaded and light shaded regions indicate 
the allowed regions combined from the two constraints taking account of the la and 2a 
statistical errors, respectively. 



